arXiv:1506.06552vl [cond-mat.mes-hall] 22Jun2015 


A Renormalization Group Study of Interacting Helical Liquid : 

Physics of Majorana Fermion 


Sujit Sarkar 

Poornaprajna Institute of Scientific Research, 

4 Sadashivanagar, 

Bangalore 5600 80, India. 

(Dated: June 23, 2015) 

Abstract 

The physics of Majorana fermion occurs at the edge and interface of many one-dimensional 
quantum many body system with gapped excitation spectrum. We present the renormalization 
group equations for strongly interacting helical liquid. We present the results of both Majorana- 
Ising topological quantum phase transition and also Berezinskii-Kosterlitz-Thouless topological 
phase transition. We show that the Majorana-Ising topological quantum phase transition for 
non-interacting case is an exactly solvable line but in presence of interaction the system has no 
exactly solvable line. We study the effect of umklapp scattering on the renormalization group flow 
diagrams and also the condition for the appearence of Majorana fermion modes. 
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I. INTRODUCTION 


The physics of Majorana fermion is the subject of intense research interest in quantum 
many body system near a deca [1-31]. Majorana fermion appeares in the topologically 
ordered state which is one of the important focous of research in condensed matter physics. 
These states are chacterized by the topological invariant quantities , Chern number and 
Winding number. Kitaev ^ and the other research group have proposed the physical 
realization of Majorana fermion as a edge state of the of one dimensional system that 
include one-dimensional superconductors, semiconductor quantum wire , proximity coupled 
to superconductor, the cold atom trapped in one-dimension. The physics and the control 
technology of Majorana fermion are the most important properties of non-abelian quantum 
computation, the quantum Hall states are the example of topological states [n|. Therefore 
the physics of Majorana fermion is almost ubiquitous in different quantum many body 
system [1-31]. 

In this research paper, we study for the existence of Majorana fermion mode in interacting 
helical liquid by using the renormalization group theory method. 

Before we proceed further we would like to state the basic properties of Helical liquid very 
briefly [9-12,16-31]. The physics of helical liquid is interesting for the following properties. 
It is generally originated from the quantum spin Hall effect in a system with or with out 
Landau levels. In the quantum spin Hall effect the left movers in the edge are connected with 
the down spin and the right movers with the up spin and the transport process is quantized. 
This physics is generally termed as a ’’Helical liquid” which describe the connection between 
the spin and momentum. It does not break time reversal invariance which occurs in chiral 
Luttinger liquid. The helical liquid has even number of time reversal pair whereas the 
spinful Luttinger liquid has odd number of pairs. The physics of helical liquid occurs in the 
surface state of a two dimensional insulator a proximity coupled semiconductor quantum 
wire subject to the spin orbit interaction. 

Here we describe the basic mathematical analysis of Helical liquid. It is well known to us 
that the low energy excitation in the one dimensional quantum many body system occurs 
near to the adjacent region of Fermi points. Therefore one can write the fermionic field 
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operator as 

= ^ ( 1 ) 

V -A<k-kF<A -A<k+kF<A 

Where A is the cut-off around the Fermi momentum (kp). We may consider the hrst term 
as a right mover [k > 0) and the second term as a left mover (/c < 0). One can write the 
fermionic held with spin a as = 'ipRa^x) + 'ipLa{x). For the low energy elementary 

excitations one can write the Hamiltonian as 


Ho 


/ dk 


( 2 ) 


The hrst term within the hrst bracket is one of the Kramer’s pair and the second term 
within the hrst bracket is the another Kramer’s pair. Therefore one of the Kramer’s pair in 
the upper edge and the other one at the lower edge of the system. The total fermionic held, 
^^{x) = ipR\ + V’Lj.j which can be expressed as a spinor -0 = {'tpLi V’Ht) • This is the simple 
mathematical picture of a helical liquid, where the spin is determined by the direction of 
the particle. The non-interacting part of the helical liquid for single edge in terms of spinor 
held is 


Hq = 'tp^wpOxa'' - ^)tp. 


( 3 ) 


The authors of Ref. 12| have introduced a model Hamiltonian to study the Majorana fermion 
in strongly interacting helical liquids. Here we describe the model very briehy of Ref. 12[ |. 
The authors have described a low-dimensional quantum many body system of topological 
insulator in the proximity of s-wave superconductor and an external magnetic held along 
the edge of this system. The additional terms in the Hamiltonian is 


5H = + B'lpLi^'ipLt + h.c. 


( 4 ) 


Where A is the proximity induce superconducting gap and B is the applied magnetic held 
along the edge of the sample. The Hamiltonian Hq is the time reversal invariant but the 
Hamiltonian 6H is not time reversal invariant. 


Now we consider the generic interaction which are time reversal, the authors of Ref. 
have considered the two particle forward and umklapp scattering as 

Hfxo = 92pJLi^i’Lii’Rt^i’Rt- 


m 

( 5 ) 
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In the umklapp scatterin g te rm, we write the umklapp term for half hlling following the 


Wu, Bernevig and Zhang 


25[ |. the umklapp term at the half hlling in a point splitted form. 


The point splited version can be described as a regularization of the theory. Therefore the 
umklapp term become 


= -gu / + a)'tpL.i\x)i^ti\x + a). 


( 6 ) 


Where a is the lattice constant. This analytical expression gives a regularized theory using 
the lattice constant a as a ultraviolet cut-off. 

If we use the hrst order Taylor series expansion of the fermionic held 


+ a) ~ adxiiR\\x). 


( 7 ) 


If we use this expansion for the umklapp scattering term then we produce the analytical 
expression for umklapp in a conventional form of the authors of Ref. 251. 


Hum = Qu'ipLi^d^^'tpLiHmdx'tpm + 


( 8 ) 


Therefore the total Hamiltonian of the system is 


H = Ho + H fug + Hum + SH. 


( 9 ) 


Now we can write the above Hamiltonian as, Hxyz = J2iHi (up to a constant) 12[ |. 


H^ = Yg - [/r + B{-iy]Si\ (10) 

Q; 

Where Jx,y = J ± A > 0 and J = vp and Jz > 0. 

Here we write down the bosonized form of this model Hamiltonian, the detail derivation is 
relegated to appendix. 


H = + K(djf) - (-)dA 

IK 71 

-cos(-\/T^0)- cosiV^TrO) + ^^cos(4\/^0) 

27r^ 


( 11 ) 


TT 71 

There is a sign mismatch between the present study and sign of Ref. [1^, please see the 
appendix for detail explanation. 


The author of Ref. 


12f| have studied the quantum phase diagram of this model Hamiltonian 
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based on the Abelian bosonization method. But here we study the same Hamiltonian 
by using the RG method for the following reasons. It is very clear from the continuum 
held theoretical study that our model Hamiltonian contains two strongly relevant and 
mutually nonlocal perturbation over the Gaussian (critical) theory. In such a situation the 
strong coupling hxed point is usually determined by the most relevant perturbation whose 
amplitude grows up according to its Gaussian scaling dimensions and it is not much affected 
by the less relevant coupling terms. However, this is not the general rule if the two operators 
exclude each other, i.e., if the held conhgurations which minimize one perturbation term do 
not minimize the other. In this case interplay between the two competing relevant operators 
can produce a novel quantum phase transition through a critical point or a critical line. 
Therefore, we would like to study the RG equation to interpret the quantum phases of the 
system. 


ANALYSIS OF RENORMALIZATION GROUP EQUATIONS 


We now study how the parameters B, A and K how under RG. through the operator 
product expansion. So the RG equations for their coefficients therefore are coupled to each 
other. We use operator product expansion to derive these RG equations which is independent 
of boundary condition 32|. In our derivation, we consider two operators, Xi = 
and X 2 = e(*“ 2 </>+*b 2 e)^ pj^ocedure, one can write these two held operators as a 

sum of fast and slow mode helds. In the fast held, the momentum range is < K < A 

and for the slow held K < Ae~'^\ where A is the momentum cut-oh, dl is the change in the 
logarithmic scale. The next step is the integration of the fast held for the operators Xi and 
X 2 , it yields a third operator at the same space time point, A 3 = e*(“i+a 2 )?i+*(bi+fe 2 ) 6 ')_ 
prefactor of A 3 can be found by the relation, A 1 A 2 ~ ^-iaia 2 +bib 2 )Hamiltonian 
consists of two operators, if we consider li and I 2 as the coefficient of the operators Ai and 
A 2 respectively. Then the RG expressions for contains the term (0102 + This 

is the procedure to derive these RG equations. 

In the RG process, one can write RG equations themselves are established in a perturbative 
expansion in coupling constant They cease to be valid beyond a certain length scale. 
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where g{l) 1 


dB 

m 

dA 

Id 

d^u 

nr 

dK 

nr 


{2-K)B + 4KguB, 

(2 - ^)a, 

(2 - 4K)g^ + 2KB^ 

4 - . 


( 12 ) 


Before we proceed further to study the detail analysis of these four RG equations. Here 
at hrst we consider the absence of umpklapp scattering, therefore the four RG equations 
reduced to three RG equations. The RG equations for the coefficients of Hamiltonian Hxyz 
are 


dA 

nr 

dB 

nr 

dK 

nr 


(2-^)A. 

(2 - K)B 

4 




(13) 


Here our main intention is to study these three equations explicitly to explain the different 
quantum phases of the system and also the Majorana-Ising quantum phase transition. After 
few steps of calculation, one can also write the above three RG equations as 


dA 

nr 

dB 

nr 

dlriK 

dl 


( 2 ->. 

(2 - K)B 
K-^A^ - KB'^ . 


(14) 


where I = ln[^] parameter, Aq is the initial value of the parameter. It is very 

clear from the above RG equations which reflect the duality in our helical liquid model 
system. The duality is the following, (p ^ 6, K K~^ and A ■(-)■ R. This duality 
in RG equations is actually from the duality, i.e., if we interchange the R and A and 
simultaneously the held 0 and 6 then the Hamitonian remains invariant of the bosonized 
Hamiltonian. 

These RG equations have trivial (A* = 0 = R*) hxed points for any arbitrary K. Apart 
from that these RG equations have also two non-trivial hxed lines, A = R and A = —R for 
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K = 1. 


It is very clear from the above RG equations, in the vicinity of K = 1, one can hnd 
the non-trivial critical points, where both of the coupling term are become the relevant 
coupling, i.e.,the coupling constant flowings off to the strong coupling limit, where each of 
them acquires a mass term. 

The coupling strengths A and B are related to the dual helds, 9 and 0 that drives the 
system to the different ground state, one is in the proximity induced superconducting gap 
and the other is in the applied magnetic held induce ferromagnetic phase Q, 3- The 
interplay between them give rise to the second order quantum phase transition at the 
intermediate couplings. We have already discussed that the K = 1, the self dual point and 
this point is also the only exactly solvable point. 

In our RG calculation, we also predict the Majorana-Ising transition . Before we proceed 
further we would like to explain the basic aspect of Majorana-Ising quantum phase transition. 


6H = + h.c. (15) 

We recast the fermionic held interms of the Majorana helds, + X 2 {x)). 

'ipRtix) = |(xi(x) -1- ix2{x)). The total Hamiltonian, H = Hq + 6H, become 

H = {iXi^d^Xi - iXi^d^Xi + imiXiXi). (16) 

Where mi _2 = A =p R (here A > 0 ). At A = R, one of the two Majorana fermion modes 
become gapless which is the signature of bulk Majorana-Ising quantum phase transition. 
We will explore this transition explicitly during the renormalization group study. 

Now we study the scaling theory for this problem. It is well known to all of us that the 
critical theory is invariant under the rescaling, the singular part of the free energy density 
satishes the following scaling relations. 

/,[A, R] = e-2'/,[e(2-V^)'A, (17) 

The scale I can be hxed from the following analytical relation, = 1 . After a few 

substituation, we arrive at the foilwing relation. 

/,[A, R] = 


7 





The phase transition of the system occurs when the coupling strength satisfy the following 
relation: 


A"(2"^)/(2-i/^)S ~ 1 (19) 

The phase boundary between this two quantum derive by using the above relation. When 
K = 1, (non-interacting case), the above phase boundary relation is A = B, which we 
have already discuss in the previous section, when we introduce the basic physics of 
Majorana-Ising topological quantum phase transition. 

Now we discuss, the effect of repulsive {K < 1) and attractive {K > 1) region of this 
phase when the umklapp interaction is present. These results are depicted in the hg.l. It 
observes from our study that the phase bounday exists for the non-interating and also for 
the interacting one. The phase bounday is a exactly solvable line for the non-interacting 
case but for the interacting case it is not exactly solvable line. 

One can understand this shift of this phase boundary from the following mathematical 
analysis. Generally one can consider the interaction between the Majorana fermions 
as Hint ~ XiXiX 2 X 2 - In the mean field level, one can do the following approximation, 
X2X2 — X2X2 >■ Thus it is clear from the mean-held analysis that one can absorab the 
effect of interaction as a redehne mass in the system (Eq.l6), which is shifting the phase 
boundary. 

Here we discuss very briehy based on this toy model about the nature of quantum phases 
during the topological quantum phase transition. 

In our toy model, we consider the simple case, i.e., the umklapp term is absent and at the 
point A = J. The complete Hamiltonian reduce to transverse Ising model. 

H, = ( 20 ) 

when B < A, the discrete Ising symmetry is spontaneously broken which yields a doubly 
degenerate ordered phase which is proximity effect induce superconducting gap state which 
form the Majorana fermion mode excitations at the edge of the system. For B > A, the 
magnetic held induce ferromagnetic state along the direction of magnetic held. 




FIG. 1: (Colour online.) This figure presents the phase boundary between the two different phases 
A and B. The color blue,magenta,red,green and yellow for K = 1,0.55,0.75,1.5,1.2 respectively. 
The blue line is the exactly solvable line. 

RESULTS AND DISCUSSIONS BASED ON RG EQUATIONS STUDY 

In fig.2, we present the RG flow diagram for A with B for K = 1. It reveals from our 
study that both the couplings ( A and B) are flowing off to the strong coupling phase. 
Here both the coupling terms are relevant, i.e., both of the coupling terms flowing off to 
the strong coupling phase. 

In hg.3, we present the RG flow diagram for A with B for K = 0.45. For this value of 
K, system is in the repulsive regime. According to the analysis of Abelian bosonization, 
the coupling term A is irrelevant but the RG flow diagram of our study shows that the 
A is non-zero but the magnitude has reduced. This is the advantage of the study of RG 
flow diagram over the Abelian bosonization study. It is also clear from this study that for 
this value of K, the magnitude of A is almost the same it’s initial value, i.e., the values 
of As are not changing during the RG flow. The proximity induce superconducting gap is 
non-zero, i.e, the Majorana fermion mode exists at the edge of topological insulator of the 
system. 

In hg.4, we present the RG flow diagram A with B for K = 0.2, i.e, the system is in the 
strongly repulsive regime compare to the K = 0.45. It reveals from our study that the 
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1 


K=1 



FIG. 2: (Colour online.) This figure shows the RG flow of A with B for non-interacting helical 
liquid, i.e, K = 1 inabsence of umklapp scattering. 


< 


K=0.45 



FIG. 3: (Colour online.) This figure shows the RG flow lines of A with B for interacting helical 
liquid in absence of umklapp scattering. Here K = 0.45 and B = 0.2. 

coupling term, A, flowing off to zero, i.e., for this case there is no existence of Majorana 
fermion mode in the system. The magnetic held induce ferromagnetic phase dieted by the 
direction of magnetic held is the only phase exist in the system. 

Therefore, it is clear from the study of hg.2, hg.3 and hg.4 that the Majorana fermion mode 
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K=0.2 



FIG. 4: (Color online.) This figure shows the RG flow lines of A with B for strongly repulsive 
regime where R = 0.2 in absence of umklapp scattering. 


< 



FIG. 5: (Colour online.) This figure shows the RG flow lines of A with B for R' = 1 in presence 
of umklapp scattering, gu = 0.2. Here B = 0.2. 

only disappear in presence of very strong repulsion , otherwise, Majorana fermion mode is 
robust for weak and intermediate values of strong repulsion. 

In fig.5 and fig. 6 , we present the results for = 0.2 and Qu = 0.4 respectively for K = 1. 
It is very clear that umklapp scattering unable to destroy the superconducting phase,i.e., 
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K= 1, g = 0.4 

’ »i| 


< 



FIG. 6: (Colour online.) This figure shows the RG flow lines of A with B for K = 1 m presence 
of umklapp scattering, Qu = 0.4. 

K= 0.45, g = 0.4 
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B 

FIG. 7: (Colour online.) This figure shows the RG flow lines of A with B for K = 0.45 and 
gu = 0.4. 

the existence of Majorana phase is always in the system. The qualitative behavious is the 
same for the both values of umklapp scattering. 

In hg.7, we present the result of the variation of A with B for K = 0.45 and Qu = 0.4. It 
is clear to us that the behaviour of the system for repulsive regime is the same in presence 
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of umklapp scattering also. The magnitude of the proximity induce superconductivity gap 
deceeasecs but it is nonzero, therefore the presence of Majorana fermion mode is still exists. 
In the present situation, our system is topological insulator in the vicinity of an s-wave 
superconductor and an external magnetic held. Therefore the appearence of the topological 
phase, i.e., the appearence of Majorana fermion mode at the edge of the topological 
insulator depends on the proximity induce gap in the system. Here we study, the length 
scale dependent of the system and show how it appears. The survival condition for the 
Majorana fermion mode at the edge of the system is the following. If the superconducting 
region has a length L, the Majorana mode survive if the condition, LA » 1, otherwise 
there is no Majorana fermion mode, i.e., the topological state of the system is absent 341. 
Therefore, it is clear from the study hg.5, £g.6 and £g.7 that the presence of umklapp 
scattering does not affect the appearence of Majorana fermion mode. 

In fig.8, we present the results of the study of A, B and K with the length scale (/). Here 
we consider three values of A = 0.1, 0.2, 0.4 as the three initial values of A but the initial 
value oi B = 0.2 and K = 1. It is clear from this length scale dependent study that the 
condition for the topological state of matter,i.e., the existence of Majorana fermion mode 
at the edge of the system. 

In fig.9, we present the results of the study of A, B and K with the length scale under RG 
process. Here we consider K = 0.45, i.e, the repulsive regime of the system. It is clear from 
our study that B is increasing very rapidtely with length scale. But for the higher values of 
As the length scale dependent is almost constant or decreases very slowly. It is clear from 
this hgure that for the smaller values of As decreasing very rapidly with length scale and 
the condition for Majorana fermion mode existence is not satisfied for the large values of A, 


34l |. and the system is in the topological 


the condition for Majorana fermion mode satisfies 
state. 

In. fig. 10, we present the variation for three different values of As with length under the 
RG process. The other parameters of the system are, K = 0.2 and B = 0.2. For this 
parameter space, the system is in the strong repulsive phase and the proximity induce 
superconducting gap goes to zero value very rapidtely. In this parameter space there is no 
evidence of Majorana fermion mode in the system. 
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FIG. 8: (Colour online.) This figure shows the variation of A,B and K with length scale. The 
curves for blue, pink and red lines are respectively for the A, K and B. A = 0.1,0.2,0.4; i? = 
0.2, K = 1, 



FIG. 9: (Colour online.) This figure shows the variation of A,B and K with length scale. The 
curves for blue, pink and red lines are respectively for the A, K and B. A = 0.1,0.2,0.4; i? = 
0.2, K = 0.45, 
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FIG. 10: (Colour online.) This figure shows the variation of A,B and K with length scale. The 
curves for blue, pink and red lines are respectively for the A, K and B. A = 0.1,0.2, 0.4; i? = 
0.2, K = 0.2, 

BEREZINSKII-KOSTERLITZ-THOULESS TRANSITION 


Berezinskii, Kosterlitz and Thouless proposed that the disordering is facilited by the 
condensation of topological defects (35|, [3^ . Generally one can consider it has as a vortex. 
It can be expressed as a two dimensional XY model. The orientation of a spin is dehne up 
to an interger multiple of 27i. One can consider the spin conhguration in which the travel 
of a closed path will see the angle rotate by 27rn where n is the topological charge. Here we 
derive and discuss the appearence of Berezinskii-Kosterlitz-Thouless (BKT) transition in 
our system. BKT transition is a topological quantum phase transition in low dimensional 
many body system. Here we derive the quantum version of it, i.e., at T = 0. 

To the best of our knowledge this is the hrst study to BKT transition in interacting helical 
liquid to hnd the physics of Majorana fermion mode. 

Before we start to discuss the appearance of BKT transition in our system, we would like 
to discuss very briefly why it is necessary to study the BKT transition. Here we study two 
different situations of our model Hamiltonian. For the hrst case, the applied magnetic held 
is absent [B = 0) and for the second case the proximity induced superconducting gap term 
is absent (A = 0). For both of these cases only one of the sine-Gordon coupling term is 
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present, therefore, there is no competition between the two mutually non local perturbation. 
Therefore one can think that there is no need to study the RG to extract the quantum 
phases and phase boundaries. But we still apply RG method for the following reason. 
Each of these Hamiltonians consist of two part, the hrst one {Hqi ) is the non-interacting 
where the 0 and 9 helds show the quadratic fluctuations and the other part of these 
Hamiltonians are the sine-Gordon coupling terms which of either 6 or (p helds. The 
sine-Gordon coupling term lock the held either 6^ or 0 in the minima of the potential well. 
Therefore the system has a competition between the quadratic part of the Hamiltonian 
and the sine-Gordon coupling term and this competition will govern the low energy physics 
of these Hamiltonians in diherent limit of the system. The diherence between the BKT 
transition and Majorana-Ising transition is that in BKT transition one transit from a 
gapped state to the gapless state (Luttinger liquid phase) but in Majorana-Ising quantum 
phase transition one transit from one gapped state to the other one. 

For the non-interacting Helical liquid system, the system shows the two sets of RG equations 
for the different limit of the parameter space. 

Here we derive the first set of BKT equations in absence of applied magnetic held {B = 0). 


dA 

m 

dK 

m 


A", 


( 21 ) 


To reduced this equations to the standard BKT equation, we do the following transforma¬ 
tions. K = 1/2K and 1 — K = —y\\- Then hnally the above RG equations become. 


dA 

m 

dy\\ 

dl 



The other set of RG equations (when A 

dB 

H 

dK 

nr 


= 0) are 
= {2-K)B 
= - K^B^ , 


( 22 ) 


(23) 


Here we do the following two transformations, K = 2 + y\\ and B —)■ B/2. Finally the 
equation reduce the standard BKT equation. 
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( 24 ) 


dB 

H 

dy\\ 

dl 

Any, BKT transition consists of three phase, the weak coupling phase, strong coupling 
phase and the intermediate coupling where the system transists from the weak coupling 
phase to strong coupling phase. The structure of these RG equations have the general 
structure for the different sets of RG equations. 

The analysis of these RG equations consists of three different phases. 

The hrst one is the weak coupling phase when > A. In this phase there is no gapped 
excitation of the system, i.e, the system is in the Luttinger liquid phase. In this phase, 
there is no evidence of Majorana fermion mode, i.e., the system is in the non-topological 
state. 

The second one is the crossover regime, the mathematical condition for this crossover regime 
is —A < y\\ < A (—A < — 1) < A). One observe the crossover from the weak coupling 

phase {y\\ = |A| ) to the strong coupling regime (|/|| = — |A|). During this phase cross over 
system transit from Luttinger liquid phase to the proximity induced superconducting gap 
phase, i.e., A 7 ^ 0. In this phase the system has the Majorana fermion mode. For this 
situation, system transit from non-topological state to the topological state. 

The third one is the deep massive phase, the RG flows flowing off to the strong coupling 
regime away from the Gaussian hxed line. 

The structure of the second BKT equations (Eq.23 ) is the same form as that of the hrst 
one (Eq. 21). Therefore the analysis of these equation are the same as that of the hrst one. 
The only diherence is that here is no evidence Majorana fermion mode in this system. The 
system only shows the ferromagnetic phase. 

The physics of BKT transition has studied in diherent condensed matter system like in one 
or two dimensional superconductivity, superhudity, Josephson junction array, melting of 
cryastaline thin him, one-dimensional metal and quantum magnets. But the explicit study 
of BKT transition for cavity QED lattice is absent in the previous literature. 
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EFFECT OF CHEMICAL POTENTIAL 


Now we consider the effect of the presence of chemical potential (/i) in the physics of 
interacting helical liqnid. At first we consider the absence of umklapp scattering [qu) term 
and the proximity indnce snpercondncting gap (A). At first we consider a transformation, 
20 —)■ 20 + 5ix, where This transformation eliminate the term dx4>{x) from 

the bosonized Hamiltonian but this transformation leads to a spatially oscillating term, i.e., 
^cos{^/%{2(|) + 5ix)) = ^cos{2^/%(|) + 5l^/%(|){x)). In the absence of gu and A, system shows 
the commensurate to incommensuarte transition. 

If 5ia >> 1, then the term ^cos{y/n{2(f> + (5ia;)) is quickly oscillating and average out to 
zero which reflect the compition between the /i and B. As a result of this compition, the 
RG flow in B has to be cut-off when 26i{l) ~ 1. To the lowest order in R, A and (5i, the RG 
flows of (5i is 


d6i 

m 


di. 


(25) 


When the all perturbations are relevant, they flows to the strong coupling phase under 
RG transformation. If the coupling B{1) reaches the strong coupling phase before 5i{l)a 
reach to one. The phase of the system is the ferromagnetic phase. This condition is 
25i{l*)a « 1 which translate into the requirement 5i(0)a << When we 

consider the opposite scenario, i.e., the B term is not able to reach the strong coupling 
phase but the A(/) term reaches the strong coupling phase and system shows the presence 
of Majorana fermion mode. The presence of has no effect on it, because the fi term is 
related with the 0 held and the proximity induce superconductivity is related with 6 term. 
Now we discuss the physics in presence of umklapp term, the system posses a oscillatory 
term ^cos(4i/^0(a;) -|- 25ii/^0(a;)). If the gu{l) term reaches the strong coupling phase 
earlier than the B{l) then the system is in the ferromagnetic phase. 


SUMMARY AND CONCLUSIONS 

We have done the renormalization group study of interacting helical liquid in presence 
of proximity induce superconductivity and applied magnetic held along the edge direction. 
We have found that Majorana fermion mode exist for this system even in the presence 
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of repulsive interaction and umklapp scattering. We have also studied the length scale 
dependence of B, A, K under the renormalization to get the condition for the existence 
of Majorana fermion explicitly. We have also studied the Berezinskii-Kosterlitz-Thouless 
transition. The condition of self duality destroy in presence of umklapp scattering. 
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Appendix 

The authors of Ref. 


12 j have not presented the detail derivation of the bosonized 


Hamiltonian. Here we present the detail derivation of the bosonized Hamiltonian during 


this derivation we follow the following references 


33 


37|. 


In the bosonization process, one can express the fermionic field of low dimensional quantum 
many body system as, 

V'R/L.t = is the Klein factor to preserve the anticommutivity of the 

fermionic field which obey the Clifford algebra (here cite reference). Here we introduce the 
two bosonic helds, 6 and 0, which are dual to each other. These two helds are related with 
the following relations, 0 = 0 k + and 6 = 6r + 6l- 


The analytical relation between the Klein factors have mentioned in Ref. 33|, l37| |. 


The non-interacting Hamiltonian for helical liquid is 

Hq = 'lpL,i\iVFd^ - + 'lpR,^^\-iVFd^ - /i)0K,t 

= VF[{d^(l)L,if + (0x0K,t)^]- ^dx4>. (26) 

Vvr 

We use the following relation during the derivation of the above Hamiltonian. 
PR/L,u — '^dx(pR/L,cr- 

Now we bosonize the forward scattering and umklapp scattering term, these two interaction 
terms are time reversal symmetry invariant terms. 


Hf — 5'20LqVL,iV’i?,tVi?,t + + (V’iidViJd) ] 

= g2PLiPR^ + ^[{PLif + [pRxf] (27) 

Now we use the bosonized version of the density operator to express the bosonized form of 
the above Hamiltonian. Finally the above Hamiltonian reduced to 


Hf — —dx(l)Lidx<t>Rt + |^[(0:r0L4.)^ + {dxfpRfY 

71 ZTl 


(28) 


Now we consider the umklapp term, we follow the convention of Wu, Bernevig and Zhang 


25| . to write the umklapp term in a point splitted form to regularize the theory where 


lattice constant use as a ultra-violet cut-off 251. 
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Hum = -Quli'ipLiHmf + 


9ii 


(27r) (27r) 

Now we derive the bosonized form of the applied magnetic held along the edge direction. 


111 


= ^VL,iVR,Ae 


_ g-2 


V 47r0^ 


(30) 


111 




1\/47v6 g—iV47r0'\ 


(31) 


= 1 


Ki = 


a^, riB^ = a" 0$ a^, Therefore one can write the Hi and k ,2 as 

\ 


V 


0 

0 

0 


0 0 

-10 0 

0-10 
0 0-1 


^ 0 0 -i 0 ^ 


Now we follow the prescription of Senechal to choose a Klein basis which simultaneously 
diagonalizes all Klein factors. Here we are faceing two Klein factors 

= VLiVRtVLiVRt- R 2 = VLiVRt- According to the notation of Senechal, we can write, 

-1 0 

I I 

0 0 0 * 

-i 0 0 0 
0 i 0 0 

The above two matrices are commute with each other, therefore one can diagonalize these 
matrices simultaneously. The eigen vectors of one matrix can be expressed as a linear 
combination of eigen vectors of other matrix from which one can construct the unitary 
matrix which diagonalize the both matrix ki and K 2 simultaneously. The diagonalize matrix. 


, «2 = 


V 


J 


we may call it as a k.^ and 

^-1 0 0 0 ^ 


K3 = U-HiU = 


V 


0 

0 

0 


-10 0 
0-10 
0 0-1 


K4 = U H 2 U = 


^ -i 0 0 0 ^ 
0 i 0 0 
0 0 i 0 
0 0 0 -i 


One can gauge choice by considering, —1 and —i from the upper left corner and forget 
the rest of the Klein space. 

Using this convention one can write the Hamiltonian Hb (Eq. 30) and H^{Eq.31) are 
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respectively. 

Hb = ^sin{y/^(f){x)). 
Ha = ^sin{y/^6{x)). 
Hgu = ^cos{y/T^(j){x)). 


v = vb + ^,K 


X — ^2 

2 'kvf * 


Explanation of signe mismatch: 

The anthor of Ref. have considered \/4^ —)■ yJAircj) + 7r/2, Att9 —)■ \/Ht9 — tt/2. 

Therefore the form of the fnll bosonized Hamiltonian is 

H = + K{djf) - 

2 K 71 

H—cos(\/T^0)- cos {V4:7r9) + -^^cos{Ay/%(p) (32) 

71 71 2n^ 
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